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£^ We review a recently developed theoretical approach to the experimental detection and 

quantification of bipartite quantum correlations between a qubit and a d dimensional 
^ system. Specifically, introducing a properly designed measure Q, the presented scheme 

^7"^ allows us to quantify general quantum correlations for arbitrary states of 2 ® d systems 

without the need to fully reconstruct them by tomographic techniques. We take in exam 
the specifics of the required experimental architecture in nuclear magnetic resonance and 
' ,. ' optical settings. Finally we extend this approach to models of open system dynamics and 

discuss possible advantages and limitations in such a context. 
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in 
od 

In recent years quantum information has been identified as an ideal ground to test 
the robustness of the fundamental principles of quantum mechanics and, at the same 
• • time, to push technology over its inherent limits. At the core of this new research 

. £^ field, a fundamental role is played by one of the key features of the quantum world, 

the superposition principle, which allows for the existence of a class of physical 
states without any classical analogue. Given a composite bipartite system, one can 
engineer quantum states of the global system in which the subsystems are not in a 
well defined state and share some amount of genuinely quantum correlations, e.g. 
entanglement. Quantum entanglement is recognized as an essential ingredient for 
most of quantum information protocols, and, therefore, a fundamental resource for 
quantum technology^. 

However, the last decade has seen a great deal of attention pointed to a subtler 
but more general kind of quantum correlations (QCs from now on). For instance, it 
has been established that separable (i.e., unentangled) mixed states can still pos- 
sess correlations of a quantum nature, which cannot be described within a classical 
probabilistic framework EES It has been conjectured that such general QCs play a 
role in performing better-than-classical algorithms for specific computational tasks 



1. Introduction 



1 
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| 5 | 6 | 7 | 8 | 9 | i ojii^ anc j th e y have been identified as a figure of merit for remote state 
preparation | 8 | 12 | 13 | and phase estimation protocols An intense experimental ac- 
tivity followed from these theoretical findingst 7 l 15 l 16 l 17 l 18 l 19 l More generally, reaching 
a full understanding of the deep nature of the QCs, finding out how we can exploit 
them, testing their accessibility and resilience under dccohcrence, and discovering 
which role they play in the study of complex and many body systems, are among the 
most exciting scientific challenges linking quantum information theory, complexity 
theory and condensed matter physics^. 

Qualitatively, general QCs are related to the disturbance induced by the mea- 
surement process on a physical systemEEll, while concepts and tools employed from 
information theory allow a quantitative evaluation of the amount of QCs in the state 
of the system. Several QCs measures have been introduced l 2 l 3 l 23 l 24 l and appealing 
operational interpretations have been associated to them^. At this stage, the nat- 
ural next step should be to establish a link between theoretical and experimental 
quantification of QCs, but, unfortunately, all the QCs measures are defined by means 
of a state-dependent optimization and are not directly associable with observable 
quantities, i.e. Hermitian operators. Since full state reconstruction is a tedious, if 
not unfeasible, procedure, it becomes desirable to find ways to evaluate QCs by 
means of a smaller number of measurements than the ones required by tomographic 
techniques. In this direction, appreciable attempts to detect nonvanishing QCs by 
observable witnesses have been realized l 15 l 16 l 18 l Anyway, reminding that almost all 
states possess QCs ^ anc l that the pivotal question around QCs is whether they 
are exploitable as resource, it seems worthy to pursue a more informative (but still 
experimentally manageable) quantitative characterization of QCs. 

In this article, we review and update the proposal two of us presented in Ref. l27L 
which bridges the gap between theoretical and experimental evaluation of QCs. 
First, QCs quantifier Q for two-qubit states is introduced, directly derived from 
the geometric discord defined in Ref. [23j which is given by a state-independent 
expression not involving any optimization procedure. Also, a generalization of Q 
to detect bipartite QCs for states of 2 (g) d dimensional systems is provided. In 
particular, as example given, we focus on the amount of QCs produced in the four- 
qubit realization of the DQC1 protocol I 5 ! 15 ! 19 !, which recently captured the interest 
of the experimental community 1 1 5 1 1 9 1 Noting that there is no need for a complete 
reconstruction of the state in order to calculate Q, we explore the possible direct 
implementation of the non-tomographic measurement scheme required to detect 
the value of Q on an unknown state. In particular, we express Q as a function of 
observable quantities (O^), i.e. the expectation values of proper Hermitian operators 
0,j. The nature and the number of such operators is obviously dependent on the 
particular setting considered. 

Two possibilities are taken in exam: first, we consider the NMR (Nuclear Mag- 
netic Resonance) setting^, in which we obtain information on the system by means 
of spin measurements. Here, our protocol allows a gain (over full state tomography) 
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which is linearly dependent on Bob's dimension. Then, we focus on the quantum 
optical setup, which implies to recast Q in terms of expectation values of projec- 
tors and swap operators. The number of measurements required is independent of 
Bob's dimension d, while the complexity of the setting, i.e., the number of opti- 
cal devices needed, increases only linearly with d. Thus, an exponential gain over 
tomography can be achieved in terms of required resources. The quantum circuits 
simulating the measurements can be designed following a well established literature 
| 29 | 30 | 3i | 32 | 33 | 34 | 35 (36[ j n g pj^ e f a m inor advantage in terms of number of measure- 
ments, the NMR implementation is by far easier to realize than the optical one 

|15|16|19[ 

The advantage provided by the correlation quantifier Q is even more striking in 
dynamical contexts. To experimentally study the evolution of QCs, many sets of 
measurements need to be performed at different instants of time, thus the number 
of measurements required to resolve the dynamics increases dramatically. The in- 
troduction of Q and of its measurement schemes then helps in reducing the number 
of necessary resources compared, for instance, to those needed for full dynamical 
tomographic reconstruction. In view of this possibility we provide a brief discussion 
on the properties of the quantum correlation quantifier Q in the context of open 
quantum systems l 37 l 38 [ Due to the prominent role that QCs promise to play for 
realistic quantum technology, in recent years strong efforts have been made to pro- 
vide a qualitative and quantitative picture on their evolution under various types of 
decoherence. The reader can find some relevant references of interest in the case of 
entanglement MMmmMEMMl, quantum discord IMMMMIM^ and other quan- 
tifiers of quantum correlations dynamics l 51 l 52 l 

Here we will focus on two different system-reservoir models for two-qubits sys- 
tems. First we consider a scenario with tunable non-Markovianity where each qubit 
is coupled with a local Lorentzian environment^, and we find that the quantifier Q 
witnesses typical non-Markovian dynamical features for quantum correlation, e.g. 
presence of oscillations due to the non-divisibility of the associated quantum dy- 
namical map^lEH The second example is instead a Markovian scenario where each 
qubit decays accordingly to a non-dissipative dynamics Recently this system 
has been invoked as a paradigmatic example showing a transition between a quan- 
tum and a classical decoherence regime, i.e. the presence of a time interval where 
quantum correlations quantified by the entropic quantum discord IHl 3 ! do not decay. 
Other quantifiers of correlations like the geometric discord ^ also witness a change 
in the decay properties of quantum correlations, despite showing a different qual- 
itative behavior. If the quantifier Q is instead employed, however, such transition 
is not evident, due to the very definition of the quantity Q which does not involve 
minimization procedures. Anyway, phenomena such as frozen discord are typical of 
states identified with a small number of parameters, e.g., Werner states, evolving 
under purpose-driven dynamics preserving such peculiarity. Obviously, the state 
reconstruction in such occurrences does not represent a serious concern, if a prior 
knowledge on the form of the states under consideration is assumed. Therefore the 
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advantages of adopting Q rather than more complex measures of QCs cannot be 
apparent in those settings. 

The paper is organized as follows. In Section [2j after a brief review on geometric 
discord Dq, we present the definition of the experimentally accessible bipartite 
QCs quantifier Q. We initially consider the paradigmatic two-qubit case, then a 
straightforward extension to 2 <Ei d dimensional systems is provided, and an explicit 
analysis carried out for the four-qubit DQC1 model ^. Section [3] describes possible 
implementations of the experimental schemes required to evaluate Q in laboratory. 
In particular, features of NMR and optical settings are explored. In Section [4] we 
focus on the open system dynamics of the QCs quantifiers Q and Dq in the two-qubit 
case, comparing their evolutions for local Markovian and non-Markovian channels. 
Finally, we summarize in Section [5] the main points of our work, with an overview 
of future developments in the field. 

2. Observable measure of QCs 
2.1. Geometric Discord 

Hereafter we focus on states p of a bipartite system AB. A quantitative assessment 
of QCs is inherently dependent on the specific measurement we are going to make on 
the system. According to a conventional choice stated by the literature of the field, 
we consider a local measurement on one of the subsystems, say Alice. The states 
left undisturbed by such a measurement are called "classical-quantum" states IHEZl, 
and form a null-measure subset f2 in the set of all density operators. Their density 
matrix takes the following form 

Pcq = ^2pt\i)(i\® PBi, (1) 

i 

where the positive coefficients pi define a probability distribution, psi are density 
operators of the subsystem B and {\i}} is an orthonormal basis for the Hilbert space 
of subsystem A. One can intuitively define the amount of QCs of a generic state in 
terms of its minimal distance from the set of classical-quantum states. Indeed, 
the geometric discord Dq, introduced for the two-qubit case in Ref. [531 is defined 
as 

D G (p) = 2 min Jp-pcqg, (2) 
pcq esi 

where the Hilbert-Schmidt norm is employed as measure of distance between states, 
i.e. ||M|| 2 = \/Tr (MAft) = rrif, and {rrii} are the eigenvalues of the matrix 

M. It should also be remarked that a normalization factor 2 is added in the definition 
|2]) , in order to obtain 1 as the maximum value for the Dq, i.e., in the case of Bell 
states. 

Geometric discord enjoys two nice theoretical interpretations. First, it quantifies 
the disturbance induced by local Von Neumann projective measurements II = 11^4 ® 
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Ib on the subsystem 

£> G (p) = 2min||p-n(p)||2. (3) 

Moreover, it can be recast as the distance of a state from itself after the action of 
a 'root-of- unity' local unitary operation on Alice U = Ua ® Ib ^ 

D G (p) = 2mm\\p-UpUi\\ 2 2 . (4) 

Recently Dq has been found out a reliable figure of merit for remote state prepa- 
ration EIHl It is appropriate to remind that D G and all the discord-like measures 
are not symmetric under subsystems swapping, i.e., performing the measurement 
on Bob rather than on Alice would lead to define another class of QCs signatures. 

Geometric discord can be easily calculated for two qubits. First, one has to write 
the state in the Bloch-Fano picturdHMl 

3 



P 



1 V - 

- 2^ Rijv 




Ujdi ®oA, (5) 



where Rij — Tr[p(ai <g> Cj)], oo = I2, <Ji (i = 1,2,3) are the Pauli matrices, x = 
{xi},y — {yi} are the Bloch column vectors associated to A,B, and ty are the 
entries of the correlation matrix t. Then one has 

D G (p) = \(\\xf +1141-4^) 

= 2(Tr[5] - fc max ), (6) 

where A: max is the largest eigenvalue of the matrix and S = j(X + T), with 
X = xx 1 \T = tt T . From now on we will imply Dq = Dc{p), assuming the same 
convention for all the other quantities. 



2.2. Observable lower bound 

The minimization remaining in the definition of Dq can be actually solved in closed 
form. Since the characteristic equation of the matrix S is a cubic with real coeffi- 
cients and roots, is easily solvable by following standard techniques ^4 Indeed, the 
eigenvalues of S can be found by solving an equation of the form 

k 3 + a k 2 +a 1 k + a 2 = 0, (7) 

where 

a = -Tr[5] 

a 1 = i(Tr[S] 2 -Tr[S 2 ]) 

a 2 = -^( ai Tr[S] + a Tr[S 2 } + Tr[5 3 ]). (8) 
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Introducing the following variables 

1 



(3oi - a ) 



1 

54 



(9a a ai - 27a 2 - 2a 3 ,) 



6 = arccos 



after a bit of algebra, one obtains 
h = i (tt[S] + 2 v /r q 3 cos 



(9) 



i(6Tr[S' 2 ] -2Tr[5] 2 ) 



= arccos 



(2Tr[S'] 3 - ^r^Tr^ 2 ] + 9Tr[5 ,3 ])i/2/(3Tr[S' 2 ] - Tr[S] 2 ) 3 



{a J = {0,27r,47r}. 



(10) 



We now have state independent expressions for the eigenvalues of S. Also, we observe 
that 9 is an arccosine function, thus its domain is < 0/3 < 7r/3. The maximum of 
cos [^p-] is therefore always reached for on = a± = 0. Hence, fc max = maxjfcj} = 
k± , and the geometric discord for an arbitrary two-qubit state p can be recast as an 
explicit function of the coefficients {pij ) 



D G = 2(Tr[5] 
2 



fci) 



2Tr[5] - v^Ti'IS 2 ] -2Tr[S'] 2 cos 



(11) 



At this point, we remind that the aim is to find an observable QCs measure, i.e. 
to quantify correlations in terms of observable quantities. The geometric discord 



in expression Eq. (11) is just a function of polynomials of the density matrix en- 



tries. Protocols for writing linear and even non-linear functionals of {pij} in terms 
of expectation values of Hcrmitian unitary operators have been extensively devel- 
oped. Furthermore quantum circuits estimating such quantities have been already 
designed For an overview of the state of the art of the field see Ref. 
Thus in principle nothing prevents us from measuring geometric discord in actual 
experimental setups. Unfortunately, in practice, the implementation of the required 
architecture seems rather challenging, hence it can be valuable to make a further 
effort and trying to define a QCs measure endowed with an even simpler and more 
accessible experimental evaluation. 



Moved by the previous considerations, we observe that in Eq. (11) one can fix 
= and define ^ 



Q=- (2Tr[S] - v/6Tr[S 2 ] - 2Tr[S] 2 ) . (12) 
It is immediate to see, from the properties of the cosine, that Q < Dq. In Fig. [T]we 
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Fig. 1. Geometric discord Dq versus Q. Sample of 10 4 randomly generated two-qubit states. The 
plotted quantities are dimensionlcss. 

compare the two quantities, showing that Q is a very tight lower bound of geometric 
discord. More important, Q is still a faithful QCs measure. Indeed, the following 
properties hold. 

• Q > 0, being zero only for classical-quantum states pcq, i.e., Q = <^=^ 
Dq — 0. To prove this, notice that the condition for vanishing Q is Tr^] 2 = 
Tr[5 2 ]. By the Cayley-Hamilton theorem, this implies Tr[S] 3 = Tv[S 3 } and 
consequently Dq — 0. 

• For pure states, Q is equal to geometric discord, as it can be easily proven. 
The Schmidt decomposition of a pure state p p of a two-qubit system reads 

Pp = V^V^K^O'ili (13) 

y=o,i 

where {cti} are the Schmidt coefficients and on = 1. Simple algebraic 
steps return 9 p = 0, thus Dq(p p ) — Q{p p ). 

In general, for two-qubit states, numerical and partial analytical evidences show 
that the following chain of inequalities holds: Dq > Q > Af 2 , where J\f is the 
negativity, a computable entanglement monotone R AH those quantities coincide 
on pure states, completing the picture presented in Ref. |63l For a more advanced 
study of the interplay between QCs and entanglement, see Ref. 19} ITT1 

2.3. Extension to 2 (g) d systems 

Here we address the problem of measuring bipartite QCs for states of 2(g) d systems, 
where subsystem A is the qubit. A generalization of geometric discord to catch 
bipartite QCs in such a case has been derived (for finite d) in Ref. [Ml Its expression 
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is the very same as Eq. ([6]). The Bloch-Fano form for the state is 




3 d'-l 



Vjh <8 Tj + 2^ 2^ UjGi ® Tj j , (14) 

i=i j=i 

where we can pick the generalized and normalized Gell-Mann matrices as basis {tj} 
of the (/-dimensional subsystem B. Obviously, {j/j} is now a d-dimensional vector 
and t is a 3 x d correlation matrix. One can notice that the matrix S = A (X + T) 
has still 3x3 dimension, thus its characteristic equation remains a cubic and we 
can repeat all the previous steps to write closed expressions for Dq and Q formally 
equivalent to Eq.( ll|12 ). The procedure can be extended to d — oo according to 



the prescription of Ref. [T3] 
2.4. DQC1 calculations 

To conclude this section, we include a simple but meaningful case study to show- 
case a comparison of QCs measures. The DQC1 algorithm, introduced in Ref. 5, 
is a non-universal quantum computing protocol estimating the trace of a rt-qubit 
unitary matrix U. Ideally suited for an NMR setting, it provides an exponential 
speed up compared to the best known classical algorithm for such a specific task 
| 6 | 7 | i5 | i9 [ rj^g sur p r igi n g feature of the DQC1 is that this enhancement in the perfor- 
mance is obtained despite a negligible amount of entanglement created during the 
computation. In particular, an ancillary qubit (Alice's subsystem) is initially in a 
state with arbitrary polarization fi, i.e., p™ = h(^-2 + M°3)> while Bob amounts to 
an n-qubit maximally mixed state, i.e., p l g = ^-la"- A four-qubit implementation 
(the ancilla A vs n = 3 qubits) has been recently investigated experimentally in 
Ref. [13 HH Specifically, the designed unitary gate is U = (a, a, b, 1, a, b, 1, 1), with 
a = — (e~' i37r / 5 ) 4 , b = (e -l3?r / 5 ) 8 . An approximation of the Jones polynomials has 
been realized by this setting ISU. The evaluation of Tr[U] runs in the following way: 
given an initial uncorrelated state p m = p™ <g) p % § , referring to the scheme of Fig. [2J 
the protocol returns the final state 

p° ut = (U H ® l B )Up m U\U H ® 1b) ] 
1 



16 



Is 






h 



(15) 



where we denote by Uh an Hadamard gate to be performed on Alice side. It is 
straightforward to see that the output state for Alice is 

,r = if '"T'Y del 



2 VfTrl^] 1 

Thus measurements of the ancilla polarization yield an estimation of the trace of 
the unitary matrix: (oi)^ = Re{Tr[[/]/8} , (a 2 ) p ont = Im {Tr[C/]/8}. 

Let us turn our attention to the correlations in the final state of the computation. 
The entanglement between the ancilla and the three-qubit register is manifestly 
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H 



II 



Fig. 2. DQC1 model with a three-qubit state in a maximally mixed state and an attached ancilla 
of purity fi. Measuring <ti,<T2 on the ancilla returns the real and imaginary part of Tr[[/], 




Fig. 3. Bipartite quantum correlations for the final state in the four-qubit DQC1 model discussed 
in the text, as measured by T> (red continuous line), Dq (blue dotted line) and Q (black dashed 
line) as functions of the initial ancilla polarization fi. All the plotted quantities are dimensionless. 



negligible and does not scale with the number of qubits involved. On the same 
hand, we can easily study QCs across this bipartition. In particular, we compare 
the behaviour of geometric QCs measures, i.e., the geometric discord Dq and the 
lower bound Q, with the entropic discord V defined in Refs. [2J [3l The former 
measures are easily calculated from Eq. ( 11|12 1 , while for the latter one we retrieve 
the approximated expression for the output states of the DQC1 model calculated 
in Ref. El 

D G = 0.0531325/i 2 
Q = 0.0402856/z 2 



V 



ho 



log 2 [l + y/l - m 2 ] - (1 - v/l-M 2 ) log 2 H, (17) 



where h 2 is the binary Shannon entropy. Surprisingly, the geometric discord of the 
DQC1 output state can be directly linked to the trace of the square of the unitary 
matrix (see Ref. for details). In Fig. [3] we study the behaviour of QCs measured 
by T>, Dq, Q by varying the purity of the ancilla in the initial state. As expected, the 
amount of bipartite QCs between the ancilla and the three-qubit register seems to 
be an efficient figure of merit for the efficiency of the protocol, as all the considered 
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QCs measures are monotonically increasing with the initial purity of the ancilla. 

3. Experimental implementation: a theoretical point of view 
3.1. Nuclear Magnetic Resonance (NMR) setting 

For this section, general references are Refs. 28, 66. In NMR experiments, quantum 
states are realized by spinorial configurations of magnetic nuclei. At room temper- 
ature, the n-qubit state of a NMR system is given by 

9 = + eAp, (18) 

where e = hiUL/2 n KBT ~ 10~ 5 is the average thermal energy. Every manipula- 
tion is implemented by varying the deviation matrix Ap, which carries the actual 
information content of the state. 

A peculiarity of the NMR setting is that there is negligible entanglement in the 
produced states. In spite of that, several quantum computational tasks have been 
satisfyingly studied and implemented by means of such a technique. It is supposed 
that QCs could be the key resources for the supraclassical performances in NMR 
environments. Indeed, the DQC1 model of computation we treated in the previous 
section was designed thinking about NMR implementation of quantum information 
processing ^. In any case, it appears a well suited ground for investigating QCs 
potentialities, and definitely this is not a coincidence. Measures of QCs such as 
geometric discord and Q are the easiest to manage at theoretical level, built by con- 
sidering the density matrix of the state in the Bloch-Fano picture. This theoretical 
framework was just introduced in the Refs. [STfl [BUI for efficiently describing the res- 
onance of magnetic nuclei under the influence of an external magnetic field. Thus, 
what seems at first sight a merely formal coincidence, underlines in fact a privileged 
interweaving between geometric quantification of QCs and NMR techniques. An 
overview of the recent studies of QCs in this setting can be found in Ref |2"81 

In the NMR context, performing global and local spin measurements is the most 
convenient method for gaining information about a state. Tomography would defi- 
nitely require the spin measurements necessary to retrieve all the Rij = Tr[(ai®Tj) p] 
coefficients, and thus the state. On the other hand, evaluating geometric discord Dq 
and the lower bound Q, by definition, does not entail to know anything about Bob's 
subsystem, i.e., we can drop the d 2 — 1 measurements related to the Bloch vector y, 
that is yj = Tr[(I 2 (8 Tj)p]. Therefore, denoting by Oj the observables to be linked 
with the QCs measures, one can write 

(Of MR ) = Tr[<7 v ® r x p\, v = 1, . . . , 3, ; A = 0, . . . , cf - 1 
D G = f[(Of MR )} 
Q = f[(Of MR )]. (19) 

Thus, it is relatively easy to quantify QCs in NMR setting l 15 l 19 l 16 l 6 H An experi- 
mental trick for further simplifying the considered procedure allows to restate global 
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spin measurements as local ones on Alice only, by just running a global operation on 
the state before the measurements^. For each particular measurement, a related 
global operation is selected. Specifically, 

€u = U C notRa((I>, n) ® n)p(U C NOTRA{(f>, n) ® Rb(4>, ")) f 

Tr \{a u ® a x )p] = Tr[(a v ® , (20) 

where we have applied local rotations R by an angle <j> about some direction n, both 
dependent on the specific a v to be evaluated, and subsequently a CNOT gate with 
Alice being the control qubit. One could maintain that when d > 2 the global spin 
measurements, i.e., estimations of the expectation values of a u ®T\ as introduced in 



Eq. (14), seem extremely intricate, and the realization of the global rotation might 
be beyond the current technological possibilities. In such a case, at least for the 
paradigmatic instance in which Bob is a n-qubit subsystem (d is even), we can pick, 
as basis {r,-} for the d-dimensional subsystem, the tensor products of Pauli matrices 

{T\} = {Id, Ci ® I<i_2,<72 ®Id_2,0"3 <g> Id-2,12 ® <Xi ® I d _ 4 , . . . , (T 3 ® . . . ® cr 3 }(21) 

reducing the detection of ijj to local spin measurements on single qubits only. In 
summary, for the NMR set up, the QCs quantification, by both geometric discord 
and Q, demands 3d 2 measurements, against the 4c? 2 — 1 required by full state recon- 
struction. Indeed, we are exempt from making local spin measurements on Bob's 
side. 



3.2. Quantum Optics 

As mentioned in the introduction, the estimation of functionals of density matrix 
elements pij has been vastly investigated by quantum optical setting. Some devices 
for the evaluation of meaningful quantities, for example the purity of the state, have 
been built having as toolbox just the very basic principles of quantum computation. 
For a broad perspective on theoretical and experimental features the reader should 
refer to Refs. M\36\ 

It has been proven that any function of the density matrix entries can be ex- 
pressed in terms of observables represented by hermitian unitary operators, and 
one can develop the experimental architecture to perform effective estimations in 
real world. It is important to stress that limits set by the actual technology could 
prevent from implementing what has been successfully designed. So, from the very 
beginning, we look for a QCs measure really accessible to experimentalists. 

Let us consider the specific case study of the implementation of Q for a two- 



qubit state. The task, as clearly expressed in Eq. (12), is to recast the quantities 
Tr[5] and Ti^S 12 ] in terms of observables. 
On can see that 

Tr[51 = ifTrLXI + TrfTl), Tr[5 2 ] = — (Tr[A 2 l + Tr[T 2 ] + 2Tr[AT]) . (22) 
4 16 
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After some algebra, one obtains 
Tr[A] = 2Tx[pi] - 1 

Tr[T] = 4(Tr[p 2 ] - Tr\p\]/2 - Tr[p|]/2) + 1 
Tr[A 2 ] = (2Tr[p A \ If 

Tr[XT] = -1 + 4Tr[p 2 ](-l + Tr[p 2 ]) + 4Tr[pi] - 4Tr[p 2 4 ] 2 + 2Tr[p|] 

+ 8Tt[ p ( Pa ® h)p(p A ® I 2 )] ~ 8Tr[p(pi ® p s )] 
Tr[T 2 ] = -32(Tr[<r 4 ] + Tr[<r 3 ]) + 3(Tr[T] 2 /2 - Tr[T] - 1/2), (23) 

where <r = p — (pa <8> l2)/2 — (I2 <8> Pb)/2. Consequently, we can write Q in terms of 
traces of multicopies of the global and marginal density matrices and their overlaps. 
In particular: 

Tr[S]=Tr[p 2 ]-Tr[pi]/2 
Tr[5 2 ] = -(-2 - 8Tr[p 4 ] + 8Tr[p 3 ] + 6Tr[p 2 ] 2 

- 2Tr[p 2 ](5 + Tr[p|]) - 2Tr[p 2 4 ] 2 + lOTr^] 

- Tr[p|] 2 + 12Tr[p|] - 6Tr[pi]Tr[p|] 

+ 4Tr[p(I 2 <g> p B )p(I 2 ® Pb)] - 24Tt[p( Pa ® ps)] 

+ 8Tr[p( Pj4 ® I 2 )p(p.4 ® I 2 )] + 8Tr[p 2 (p A ® p B )]). (24) 



Therefore, Q can be recast as a functional of polynomials (of up to the fourth order) 
of the density matrix elements, specifically traces of matrix powers and overlaps. We 
identify nine independent terms in the expressions of Eq. ( 24 ) . Inspired by the histor- 
ical lesson of nontomographic entanglement detectionl 68 ! 69 !, we associate to them the 
expectation values of the operators {OP PT }i =1 - Explicitly, the Of PT are swap/shift 
operators V k acting on k (k < 4) copies of the global and/or marginal density ma- 
trices and related overlaps. For a density matrix p it holds that Tr[p fc ] = Tr[V fe p® fc ] 

|29|30|31|32|35|36^ where yk jg the shift operatorj V k \^ 2 ■■■^k) = IV'fc^l ■ ■ ■ ^fe-l)- 

Also, for two unknown states p%, P2, it has been proven that Tr[y 2 pi<g)p 2 ] = Tr[pip 2 ] 
| 32 | 33 | 34 [ ]yj ore generally, we have Tr[pip 2 . . . pk] = Tv[V k pi® p 2 ® ■ . - ®pk]- We briefly 
present a proof of the last statement, see also Ref. [7(2 for a more elegant treatment. 
One can see that Tr[pip 2 ■ ■ ■ pk] = Tr[F fc pi<g>p 2 <g). . .®pk] by expanding the left-hand 
term in the equation as follows: 



Tr[p 1 p 2 ...p fe] = x pr J1 p? J2 ---Pk-i 3k _ 1 Pk k 3 jiw:...sir^ 

= E Pi^P2 l3 ---PktPkn- (25) 



Denoting by {\i)} a Hilbert space basis, px...k — ®iPi> an d building the shift operator 
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|0) 



Pk 



— H 




H 































Fig. 4. Circuit estimating Tr[OiPi (8> pi ® ■ ■ ■ <8> Pk] — v. A Hadamard gate H is applied to the 
ancillary qubit, followed by a controlled-Oi gate acting on the overlap of states and then another 
Hadamard one. Then, a measurement in the computational basis returns the visibility v, i.e., the 
expectation value of the operator 0.p pT . 



as chain of swaps V 2 = J2i A i B j A j B luisX-M^sli one has 

V k pi (g) p 2 ® Pk= Pl.'.'.kji...j k Niifc • ■ ■ *fc-i) Uik ■ ■ ■ jk-i \ ■ ■ ■ 

\h . . .ik-ijk){ji ■ ■ -jk-ihlh ■ ■ ■ ik-iik){ji ■ --jk-dkl, 

(26) 

while the right-hand term of the initial relation takes the form 
Tr [V k Pl ® p 2 ® . . . ® p k ] = PtA.-.iX 

ii...i k 

= E pi\A 2 i 3 --pktpi k i 1 i (27) 

ii— ik 

thus the assertion is proven. For example, (OP PT ) = Tr[p 4 ] = Tr[V 4 /9® 4 ], and so 
forth for the other terms. All the quantum circuits to be implemented for esti- 
mating (OP PT ) (by inserting an ancillary qubit) have the same architecture, which 
is depicted in Fig(4j A Mach-Zender interferometer is modified by a controlled-C^ 
gate. At the end of the routine, the visibility v related to the interference fringes 
yields the expectation value of the operator Of PT on the general dummy overlap 
px ® p2 ® • • • <& pk, as one can write Tr[Op PT pi ® p 2 ® . . . ® pk] = v. Hence, to 
quantify the degree of QCs for an arbitrary two-qubit state p, this method requires 
9 independent measurements instead of the 15 necessary for tomography^. 

Anyway we can indeed do better than this. In particular, also looking at Ref. l6"2l 
we may appreciably reduce the number of measurements and the complexity of the 
setting. We observe that V 2 = Vij = Y^ij lu)C?'*l = |(^4+Sfe ^k^Ck) (note that the 
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V defined in Ref. 1621 is twice ours). Then, let us introduce the following quantities 



Cl 


= T±[(PX iAa ( 






C2 


= n(p Al A 2 < 


*>I Bl B 2 )(p® 2 )} 




C3 


= Tr[(I AlA2 $ 






C4 


= Tr[(P AlA4 ( 


*> ^A 2 A 3 ® P Bl B 2 < 


^B 3 B 4 )(^ 4 )] 


c 5 


= Tr[(P^ lA4 ( 


E) Ia 2 a 3 ® f£ lB2 € 


^B 3B4 )(P 04 )] 


C6 


= Tr[(PI lA4 ( 


S ^A 2 A 3 ® ^B 2 ( 




C7 


= Tr[(I AlAi $ 




^b 3 bJ(p®% (28) 




Vij) is the projector on the antisymmetric subspace for a two- 



where P y = |(1 

qubit state. Evaluating Cj is equivalent to measuring the quantities 

di =Tr[(K4 lA2 ®y BlB2 )(p 82 )] 
d 2 =Tr[{l AlA2 ®V BlB2 )(p® 2 )] 
d 3 = Tr[(I AlAi ® <g> T/ SlB2 ® ^B 3 B 4 )(p® 4 )] 

d± = Tr[{V AlA4 ® Vk 2 A 3 ® Vb iB2 ® V^ 3B4 )(p^ 4 )]. (29) 
More important, one can see that 



Tr[S] = 4ci - 2c 2 - c 3 + ^ = di - ^ 2 
Tr^ 2 ] = 16c 4 + 8(c 7 - c 5 - 2c 6 ) + c 2 + 4c^ - c 3 - 2c 2 + - 



^4-^3 + 7^2- ( 30 ) 



1 

4 

Therefore, Q could be detected by evaluating seven projective or even just four 
swap measurements: 

Q = 9l(0? PT )} 
{0° PT } = {a} or {di}. (31) 

We remark that geometric discord would require a rather more complex expression 
in terms of overlaps or alternatively measurements over six copies of the state 
entailing by far a harder implementation. 

Now, let us have a look at the extension to the 2 ® d case. We can arguably 
say that the very same expressions hold, at least at formal level. Clearly, we have 
to generalize the swap and the projectors to arbitrary finite dimension. A state of 
a d-dimensional system reads p = h(^d + J2i x i T i)i implying Tr[p 2 ] = k(l + \x\ 2 ). 
Thus, in the most general fashion one obtains 

V = Jo^+XVjBTi), (32) 



and consequently 



P~ = ^((d-i)I d 2-J2n®n), (33) 
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where the TjS reduce to Pauli matrices in d = 2. The optical implementation of pro- 
jectors on PsiBj, i-e. multiqubit projectors, is more complicated than the two-qubit 
case, see e.g. Ref. |TJ However, the method presented in Ref. [7T] can be extended to 
arbitrary dimensions without terribly increasing the complexity of the experimen- 
tal setting. More precisely, the necessary number of optical devices should increase 
polynomially with d, not exponentially. See also Ref. [72] for more a detailed analy- 
sis. One can appreciate that the number of measurements is independent of Bob's 
dimension. Anyway, it could be pointed out that swaps and projectors on large di- 
mensional systems seem of hard implementation. However, such a question can be 
easily overcome in two ways. First, as already done for the NMR setting, we can 
restrict to the case of even d and pick for Bob's subsystem the basis considered in 



Eq. (21 1, such that V,P~ can be easily rewritten in terms of two-qubit operators. 
In such a case, the number of measurements would increase linearly with d not 
compromising the scalability of the protocol. Alternatively, we can directly rewrite 
d-dimensional swaps/projectors as matrix products of their two-qubit versions. We 
leave for future investigations an extensive treatment of this issue. 

4. Quantum correlation dynamics in open systems 
4.1. Overview 

In this section we will have a look at the properties of QCs, and in particular of 
the quantifier Q, in the context of dynamical open quantum systems. The general 
aim of the study of open quantum systems is to provide tools to characterize the 
dynamics of the state ps{t) of a system of interest S when it interacts with an 
environment E l 37 l 38 [ One of the possible ways to solve this problem is to derive 
an equation of motion for ps(t), which is usually obtained starting from the Von 
Neumann equation for the global state of system and the environment, then tracing 
out the degrees of freedom of the environment. The result is an operator equation 
for the state of the system only, whose structure depends however on the system- 
environment interaction and on all the physical parameters involved, e.g. frequencies 
of the environmental degrees of freedom or coupling strengths. 

The dynamics of correlations in open system contexts, e.g. entanglement or 
more general QCs, has been analyzed deeply and broadly in the last years 
|39|40|4i|42|43|44|47|48|45|49|50i A key effect of the interaction with an environment is, 

for instance, the loss of bipartite correlations for long times, with the consequence 
of a final equilibrium state not presenting typical signatures of quantumness. While 
quantum features are lost during the whole evolution, it is not uncommon to observe 
a partial restoration of such properties for certain time intervals of the dynamics. 
These non-Markovian effects arc interpreted as a consequence of memory effects of 
the environment, capable of restoring quantum coherence or correlations for short 
periods of time. Such a behavior usually is manifested in the very early stages of 
the dynamics (non-Markovian time scales) at least for weakly coupled systems, e.g. 
optical systems, and therefore it is experimentally challenging to observe. A recent 
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successful attempt is reported in Ref. [75] where the authors make use of the measure 
introduced in Ref. 1541 Another interesting result appears in Ref. 1561 where the au- 
thors show how QCs measured by quantum discord are frozen in the early stages of 
the dynamics. Such a transition from a classical to a quantum decoherence appears 
qualitatively different if other measures of quantum correlations are employed, e.g. 



geometric discord^. As we will review in Section 4.3 geometric discord does not re- 
main constant at any time, but a change of dynamical behavior is anyway observed 
at the transition time evaluated in Ref. 1561 

The main point of interest of this section is to establish whether the approximate 
value of geometric discord captured by the quantifier Q is faithful during the open 
system dynamics of the system of interest. In the two examples we consider in the 
next section we will see that, since the state changes in time according to the open 
system dynamics, the difference between Dq and Q will be in general function of 
time. Differences in the dynamics of Dq and Q are clearly due to the peculiar 



evolution of the parameter 9 defined in Eq. ( 10 ) 



4.2. Independent environments: non-Markovian case 

The first open system model we introduce is a system made of two identical non- 
interacting qubits, each coupled with its own bosonic reservoir at zero temperature. 
The total Hamiltonian H is then written as a sum of two terms, Hi and H2 having 
the form 

H t = wo^VM +5> fc6 «t 6 M +5>(4 i >6« +a«6£>t) (34) 
fc fc 

where ujq is the frequency of the qubits, the u>i- are the frequencies of the envi- 
ronmental oscillators, the gk quantify the coupling of the oscillators with the fc-th 
mode of the respective environment, Oj! and are the Pauli raising and lowering 
operators for the i-th oscillator, and, finally, the bj^ and b^ are the annihilation 
and creation operators for the mode k of the i-th reservoir. From the previous 
definition it should be clear we are under the assumption that the environments 
have the same properties and are coupled with the same strength to the respective 
qubits. The properties of the reservoirs at zero temperature can be all condensed in 
the knowledge of their spectral function J(u>). Here we choose a Lorentian shaped 
spectral distribution of the form 

JH = 27 ( Wo + A2 (35) 

where 70 is an overall coupling constant, A is the width of the spectral function and 
ujq is the center of the spectrum, which is resonant to the free qubit frequencies. 

Since the total Hamiltonian is given by the sum of two commuting terms, for any 
initial state the dynamics of the system can be solved directly from the knowledge 
of the solution for a single qubit in its environment as proven in Ref. 1531 Given the 
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density operator p(0) at the initial time to = 0, its density matrix elements in the 
canonical basis {|1) = |11>, |2> = |10), |3) = |01), |4) = |00)} at time t are given by 





= Pn(0)P t 2 , 




P22{t) 


= p22(0)P t +Pll{0)P t (l-Pt) 




= p 33 (0)P t + Pii(0)P t (l- 


Pt), 


pa{t) 


= 1 - 0Ll(f) - 022(f) - /933(f), 


P12(*) 


= Pi 2 (0)P t 3/2 , 




Pis(f) 


= Pl3(0)P t 3/2 , 


Pli(t) 


= Pu(0)Pu 




023 (*) 


= 02 3 (O)P t , 


P2i(t) 


= P t 1/2 [ P2 4(0) + p 13 (0)(l- 


-Pt)} 


P34(t) 


= P t 1/2 [ P34 (0) + p 12 (0)(l-P t 



where we defined 



Pt 



cos Ai 



2A 



sin At 



(36) 



(37) 



with A = y/2joX — A 2 /2. Equation (37) contains all the information on the dy- 
namics. Since we have the expression for the state we can now apply the results 
of the previous sections and study the dynamics of the geometric discord Dq ^ 
and quantifier Q ^ for any initial state. Here we consider the case of initial Werner 
states, having the form 



Pw(r)=r\+)(+\ + 



1-r 1 

i I4 = T 



/ 1 - r \ 
1 + r 2r 
2r 1 + r 
V 1-r/ 



(38) 



where | +) = ( 1 2) + 1 3) )/ \f2 is a maximally entangled Bell state. During the dynamics 
described by Eqs. ( [36| and ( [37] ) the evolution brings the system through states for 
which the geometric discord Dq and the quantifier Q assume different values. Just 



as an example we show in Fig. 5(a) the dynamics for an initial Werner state with 
r = 1/2 as a function of t. The environment is characterized by 70 = 1, and two 
different values for the width, A = 1 (Markovian case) and A = 0.1 (non-Markovian 
case). In both situations, Dq and Q show the same qualitative dynamics (decay 
and oscillations) . A small difference in the actual amount of correlations quantified 
appears in the first stage of the evolution. For completeness we show in Fig. [5] the 
maximum value of this difference, max t> o[-D(3(£, r ) — Q{^i r )\ as a function of the 
Werner parameter < r < 1. The result clearly states that as the parameter r 
grows, the quantifier Q is a less tight bound. 



4.3. Depolarizing Markovian channel 

The second example we consider here is the case of two qubits evolving accordingly 
to two identical independent depolarizing channels The equation of motion for 
the single qubit channel has the form 

-Pi^vfpi^-Pi) (39) 
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r 



(a) (b) 



Fig. 5. (a) Correlation quantifier Q and geometric discord Dq for an initial two qubit Werner 
state with r = 3/4 as a function of time t. Environment parameters 70 = 1, A = 1 (blue solid, 
Q and red dashed Dq lines) and A = 0.1 (black dotted, Q, and green dotted-dashed, Dq, lines), 
(b) Maximum difference between geometric discord Dq(1>) and quantifier Q(t) for different initial 
Werner states (parameter r is varied). 



where i = A, B is the qubit index, while c- , af^ and cr| are three Pauli matrices 
for the i-th qubit. For more details the reader can refer to 1561 and references therein. 
We consider here only the phase flip channel case, i.e. j = 3, and initial maximally 
mixed marginals states of the form 

PAB = \ (lAB + E C i(°)"i i °i l ) ( 4 °) 
^ i=l ' 

with < |cj| < 1. The time evolution of the state is given by 

PAB® = A+|*+)(*+| +A+|$+)($+| +A $ |$-)($-| +A^|*-><*-| (41) 

where |*+) = (|00) ± |11))/V2 and |$+) = (|01) ± |10))/\/2 are Bell states and 

A± = [l± Cl (i)Tc 2 (<) + c 3 (t)]/4 
A± = [l± Cl (t)±c 2 (t)-c 3 (t)]/4 

Finally we have c\(t) = ci(0) exp(— 2^t), c^it) = 02(0) exp(— 2jt) and 03(f) = 03(0). 



In Fig. 6(a) we plot the dynamics of geometric discord Dq and the quantifier 
Q for 7 = 1, ci(0) = 1 and c 3 (0) = — c 2 (0) = 0.6. As one could expect, the 
sudden transition is witnessed only by the geometric discord, while the evolution 
of the quantifier does capture this dynamical feature. From Fig. [6] we show the 
maximum difference between geometric discord and Q. This maximum happens to 
be in correspondence with the transition point. 

The results of this section show that, essentially, also in dynamical contexts (at 
least in these two examples shown) the correlation quantifier Q behaves as a good 
tight lower bound to the actual geometric discord and therefore its use is encouraged 
also in the case of open system dynamics thanks to its improved accessibility. Some 
distinctive dynamical features are however not captured, and the origin of this has to 
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Q,D G 




Fig. 6. (a) Time evolution of the correlation quantifier Q (blue solid line) and geometric discord 
Dq (red dashed line) for an initial two qubit state of the form l |40[ l with ci (0) = 1, C3 (0) = — C2 (0) = 
0.6 and 70 = 1. (b) Maximum difference between geometric discord DQ(t) and quantifier Q(t) for 
different initial states (parameters £3(0) = — C2(0) are varied). 



be linked to the fact that, unlike Dq, Q is defined without a minimization procedure 
over projective measurements. 



5. Conclusions 

We discussed a proposal for quantifying bipartite QCs in states of 2(g) d dimensional 
composite systems at both theoretical and experimental level, building upon the 
results of Ref. 27. An analytical study of geometric discord led to define a tight 
meaningful lower bound Q for two-qubit states, successively extended to 2 ® d 
dimensions. From a theoretical point of view, Q is faithful and easy to calculate, 
being characterized by a state-independent expression. On the experimental side, 
Q is the friendliest measure we have so far, as its evaluation in laboratory appears 
to be in the reach of current technological levels and flexible enough to allow direct 
implementations in both NMR and optical settings. Indeed, Q is just a function of 
polynomials of the density matrix entries, thus we can set a number of operators 
{Oi} whose expectation values yield the value of Q, i.e., the amount of QCs of the 
state. On this purpose, we have presented the accessible methods and the necessary 
and sufficient resources to detect the amount of QCs of an unknown 2 eg) d state for 
both the setups. In particular, in the NMR implementation, deeply entwined with 
the geometric perspective for QCs as we discussed in this work, the detection of 
Q demands a number of local spin measurements smaller than tomography, with 
such a gain growing linearly with Bob's dimension. On the other hand, considering 
optical devices, we identified the circuits to be implemented for the measurements 
of {(Oi)}: swap or projective measurements over up to four copies of the state are 
required to detect QCs. 

It is known that a single observable is sufficient to detect nonvanishing QCs 
|!5|i6|i7|23[ Q n ^- ne othgj- hand, in general, we need a set of at least four independent 
measurements to quantify the amount of QCs for states of 2 <g> d systems. The 
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dimensionality of the Hilbert space of the system (in particular, the knowledge that 
subsystem A is a qubit) is the only information we need a priori from the state, 
and this can be retrieved by a limited amount of supplemental resources ^3 

The advantage of using Q as a quantifier of QCs is even more evident when 
we are interested in dynamical contexts where the effort to reconstruct the state 
and therefore the correlations using standard tomographic approaches becomes even 
more demanding. Being a tight lower bound to geometric discord, Q captures essen- 
tially the same features under decoherence dynamics, e.g. oscillations, sudden death 
and revivals. However attention must be payed in certain cases, e.g. sudden transi- 
tions in the dynamics of correlations, where instead Q may manifest qualitatively 
different, typically smoother features. 

Reaching tangible advances on the theoretical and experimental characterization 
of general QCs is relevant for quantum information theory, quantum foundations 
and the study of complex and many-body systems. We hope to have achieved a 
useful step forward in this scientific effort by providing new insights on QCs quan- 
tification. Nature is pervaded by evident and hidden manifestations of quantumness, 
and it seems legit to conjecture that some of its fundamental regulative processes 
are ruled and can be understood by electing QCs as primary tool of investigation. 
For example, we could envisage that the QC dynamics in a bipartite system S 
is somewhat entwined with decoherence produced by the interaction of the sys- 
tem with an environment E. Furthermore, even if it is still not conclusively known 
whether QCs in absence of entanglement can provide some actual appreciable help 
in developing better-than-classical protocols for any quantum information task, it 
is undoubtedly worthy to explore such a possibility. Experimentally sound recipes 
for measuring QCs for large dimensional systems, and serious engagement in study- 
ing their dynamics in the context of open quantum system, are pivotal to assess 
the role of QCs and their usefulness for the performance of such practical tasks in 
realistic conditions and, at the same time, to shed light on foundational questions 
of the broadest scientific prominence, such as the ultimate quantum picture of the 
measurement process. 
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